Genomic insights into the secondary aquatic transition of penguins.

Penguins lost the ability to fly more than 60 million years ago, subsequently evolving a hyper-specialized marine body plan. Within the framework of a genome-scale, fossil-inclusive phylogeny, we identify key geological events that shaped penguin diversification and genomic signatures consistent with widespread refugia/recolonization during major climate oscillations. We further identify a suite of genes potentially underpinning adaptations related to thermoregulation, oxygenation, diving, vision, diet, immunity and body size, which might have facilitated their remarkable secondary transition to an aquatic ecology. Our analyses indicate that penguins and their sister group (Procellariiformes) have the lowest evolutionary rates yet detected in birds. Together, these findings help improve our understanding of how penguins have transitioned to the marine environment, successfully colonizing some of the most extreme environments on Earth.

Adaptation and Cryptic Pseudogenization in Penguin Toll-Like Receptors.

Penguins (Sphenisciformes) are an iconic order of flightless, diving seabirds distributed across a large latitudinal range in the Southern Hemisphere. The extensive area over which penguins are endemic is likely to have fostered variation in pathogen pressure, which in turn will have imposed differential selective pressures on the penguin immune system. At the front line of pathogen detection and response, the Toll-like receptors (TLRs) provide insight into host evolution in the face of microbial challenge. TLRs respond to conserved pathogen-associated molecular patterns and are frequently found to be under positive selection, despite retaining specificity for defined agonist classes. We undertook a comparative immunogenetics analysis of TLRs for all penguin species and found evidence of adaptive evolution that was largely restricted to the cell surface-expressed TLRs, with evidence of positive selection at, or near, key agonist-binding sites in TLR1B, TLR4, and TLR5. Intriguingly, TLR15, which is activated by fungal products, appeared to have been pseudogenized multiple times in the Eudyptes spp., but a full-length form was present as a rare haplotype at the population level. However, in vitro analysis revealed that even the full-length form of Eudyptes TLR15 was nonfunctional, indicating an ancestral cryptic pseudogenization prior to its eventual disruption multiple times in the Eudyptes lineage. This unusual pseudogenization event could provide an insight into immune adaptation to fungal pathogens such as Aspergillus, which is responsible for significant mortality in wild and captive bird populations.

Independent evolution of intermediate bill widths in a seabird clade.

Interspecific introgression can occur between species that evolve rapidly within an adaptive radiation. Pachyptila petrels differ in bill size and are characterised by incomplete reproductive isolation, leading to interspecific gene flow. Salvin's prion (Pachyptila salvini), whose bill width is intermediate between broad-billed (P. vittata) and Antarctic (P. desolata) prions, evolved through homoploid hybrid speciation. MacGillivray's prion (P. macgillivrayi), known from a single population on St Paul (Indian Ocean), has a bill width intermediate between salvini and vittata and could also be the product of interspecies introgression or hybrid speciation. Recently, another prion population phenotypically similar to macgillivrayi was discovered on Gough (Atlantic Ocean), where it breeds 3 months later than vittata. The similarity in bill width between the medium-billed birds on Gough and macgillivrayi suggest that they could be closely related. In this study, we used genetic and morphological data to infer the phylogenetic position and evolutionary history of P. macgillivrayi and the Gough medium-billed prion relative other Pachyptila taxa, to determine whether species with medium bill widths evolved through common ancestry or convergence. We found that Gough medium-billed prions belong to the same evolutionary lineage as macgillivrayi, representing a new population of MacGillivray's prion that originated through a colonisation event from St Paul. We show that macgillivrayi's medium bill width evolved through divergence (genetic drift) and independently from that of salvini, which evolved through hybridisation (gene flow). This represents the independent convergence towards a similarly medium-billed phenotype. The newly discovered MacGillivray's prion population on Gough is of utmost conservation relevance, as the relict macgillivrayi population in the Indian Ocean is very small.

Vertical distribution of prokaryotes communities and predicted metabolic pathways in New Zealand wetlands, and potential for environmental DNA indicators of wetland condition.

Globally, wetlands are in decline due to anthropogenic modification and climate change. Knowledge about the spatial distribution of biodiversity and biological processes within wetlands provides essential baseline data for predicting and mitigating the effects of present and future environmental change on these critical ecosystems. To explore the potential for environmental DNA (eDNA) to provide such insights, we used 16S rRNA metabarcoding to characterise prokaryote communities and predict the distribution of prokaryote metabolic pathways in peats and sediments up to 4m below the surface across seven New Zealand wetlands. Our results reveal distinct vertical structuring of prokaryote communities and metabolic pathways in these wetlands. We also find evidence for differences in the relative abundance of certain metabolic pathways that may correspond to the degree of anthropogenic modification the wetlands have experienced. These patterns, specifically those for pathways related to aerobic respiration and the carbon cycle, can be explained predominantly by the expected effects of wetland drainage. Our study demonstrates that eDNA has the potential to be an important new tool for the assessment and monitoring of wetland health.

Latitudinal, sex and inter-specific differences in mercury and other trace metal concentrations in Adélie and Emperor penguins in the Ross Sea, Antarctica.

We sought to determine mercury (Hg) and other trace metal concentrations in Adélie (Pygoscelis adeliae) and emperor penguin (Aptenodytes forsteri) breast feathers from the Ross Sea, Antarctica, and relate those concentrations to the trophic position and the habitats in which each of these species forage. Adélie penguin feathers from the southern Ross Sea colonies were higher in Hg than those sampled further north in the Ross Sea, potentially due to greater exposure to local sources, such as volcanism. Female Adélie penguins had lower feather total Hg concentrations than males. This may reflect female penguin's capacity to eliminate Hg through the egg development and laying process, or the larger and/or older prey items that male birds can consume, reflected by their higher trophic position. Emperor penguins have higher Hg concentrations than Adélie penguins which is also partially explained by Adélie penguins feeding at lower trophic levels than emperor penguins.

Receding ice drove parallel expansions in Southern Ocean penguins.

Climate shifts are key drivers of ecosystem change. Despite the critical importance of Antarctica and the Southern Ocean for global climate, the extent of climate-driven ecological change in this region remains controversial. In particular, the biological effects of changing sea ice conditions are poorly understood. We hypothesize that rapid postglacial reductions in sea ice drove biological shifts across multiple widespread Southern Ocean species. We test for demographic shifts driven by climate events over recent millennia by analyzing population genomic datasets spanning 3 penguin genera (Eudyptes, Pygoscelis, and Aptenodytes). Demographic analyses for multiple species (macaroni/royal, eastern rockhopper, Adélie, gentoo, king, and emperor) currently inhabiting southern coastlines affected by heavy sea ice conditions during the Last Glacial Maximum (LGM) yielded genetic signatures of near-simultaneous population expansions associated with postglacial warming. Populations of the ice-adapted emperor penguin are inferred to have expanded slightly earlier than those of species requiring ice-free terrain. These concerted high-latitude expansion events contrast with relatively stable or declining demographic histories inferred for 4 penguin species (northern rockhopper, western rockhopper, Fiordland crested, and Snares crested) that apparently persisted throughout the LGM in ice-free habitats. Limited genetic structure detected in all ice-affected species across the vast Southern Ocean may reflect both rapid postglacial colonization of subantarctic and Antarctic shores, in addition to recent genetic exchange among populations. Together, these analyses highlight dramatic, ecosystem-wide responses to past Southern Ocean climate change and suggest potential for further shifts as warming continues.

High-coverage genomes to elucidate the evolution of penguins.

BACKGROUND: Penguins (Sphenisciformes) are a remarkable order of flightless wing-propelled diving seabirds distributed widely across the southern hemisphere. They share a volant common ancestor with Procellariiformes close to the Cretaceous-Paleogene boundary (66 million years ago) and subsequently lost the ability to fly but enhanced their diving capabilities. With ∼20 species among 6 genera, penguins range from the tropical Galápagos Islands to the oceanic temperate forests of New Zealand, the rocky coastlines of the sub-Antarctic islands, and the sea ice around Antarctica. To inhabit such diverse and extreme environments, penguins evolved many physiological and morphological adaptations. However, they are also highly sensitive to climate change. Therefore, penguins provide an exciting target system for understanding the evolutionary processes of speciation, adaptation, and demography. Genomic data are an emerging resource for addressing questions about such processes. RESULTS: Here we present a novel dataset of 19 high-coverage genomes that, together with 2 previously published genomes, encompass all extant penguin species. We also present a well-supported phylogeny to clarify the relationships among penguins. In contrast to recent studies, our results demonstrate that the genus Aptenodytes is basal and sister to all other extant penguin genera, providing intriguing new insights into the adaptation of penguins to Antarctica. As such, our dataset provides a novel resource for understanding the evolutionary history of penguins as a clade, as well as the fine-scale relationships of individual penguin lineages. Against this background, we introduce a major consortium of international scientists dedicated to studying these genomes. Moreover, we highlight emerging issues regarding ensuring legal and respectful indigenous consultation, particularly for genomic data originating from New Zealand Taonga species. CONCLUSIONS: We believe that our dataset and project will be important for understanding evolution, increasing cultural heritage and guiding the conservation of this iconic southern hemisphere species assemblage.

More than the eye can see: Genomic insights into the drivers of genetic differentiation in Royal/Macaroni penguins across the Southern Ocean.

The study of systematics in wide-ranging seabirds can be challenging due to the vast geographic scales involved, as well as the possible discordance between molecular, morphological and behavioral data. In the Southern Ocean, macaroni penguins (Eudyptes chrysolophus) are distributed over a circumpolar range including populations in Antarctic and sub-Antarctic areas. Macquarie Island, in its relative isolation, is home to a closely related endemic taxon - the royal penguin (Eudyptes schlegeli), which is distinguishable from E. chrysolophus mainly by facial coloration. Although these sister taxa are widely accepted as representing distinct species based on morphological grounds, the extent of their genome-wide differentiation remains uncertain. In this study, we use genome-wide Single Nucleotide Polymorphisms to test genetic differentiation between these geographically isolated taxa and evaluate the main drivers of population structure among breeding colonies of macaroni/royal penguins. Genetic similarity observed between macaroni and royal penguins suggests they constitute a single evolutionary unit. Nevertheless, royal penguins exhibited a tendency to cluster only with macaroni individuals from Kerguelen Island, suggesting that dispersal occurs mainly between these neighboring colonies. A stepping stone model of differentiation of macaroni/royal populations was further supported by a strong pattern of isolation by distance detected across its whole distribution range, possibly driven by large geographic distances between colonies as well as natal philopatry. However, we also detected intraspecific genomic differentiation between Antarctic and sub-Antarctic populations of macaroni penguins, highlighting the role of environmental factors together with geographic distance in the processes of genetic differentiation between Antarctic and sub-Antarctic waters.

Additive Traits Lead to Feeding Advantage and Reproductive Isolation, Promoting Homoploid Hybrid Speciation.

Speciation through homoploid hybridization (HHS) is considered extremely rare in animals. This is mainly because the establishment of reproductive isolation as a product of hybridization is uncommon. Additionally, many traits are underpinned by polygeny and/or incomplete dominance, where the hybrid phenotype is an additive blend of parental characteristics. Phenotypically intermediate hybrids are usually at a fitness disadvantage compared with parental species and tend to vanish through backcrossing with parental population(s). It is therefore unknown whether the additive nature of hybrid traits in itself could lead successfully to HHS. Using a multi-marker genetic data set and a meta-analysis of diet and morphology, we investigated a potential case of HHS in the prions (Pachyptila spp.), seabirds distinguished by their bills, prey choice, and timing of breeding. Using approximate Bayesian computation, we show that the medium-billed Salvin's prion (Pachyptila salvini) could be a hybrid between the narrow-billed Antarctic prion (Pachyptila desolata) and broad-billed prion (Pachyptila vittata). Remarkably, P. salvini's intermediate bill width has given it a feeding advantage with respect to the other Pachyptila species, allowing it to consume a broader range of prey, potentially increasing its fitness. Available metadata showed that P. salvini is also intermediate in breeding phenology and, with no overlap in breeding times, it is effectively reproductively isolated from either parental species through allochrony. These results provide evidence for a case of HHS in nature, and show for the first time that additivity of divergent parental traits alone can lead directly to increased hybrid fitness and reproductive isolation.

Mitogenomes Uncover Extinct Penguin Taxa and Reveal Island Formation as a Key Driver of Speciation.

The emergence of islands has been linked to spectacular radiations of diverse organisms. Although penguins spend much of their lives at sea, they rely on land for nesting, and a high proportion of extant species are endemic to geologically young islands. Islands may thus have been crucial to the evolutionary diversification of penguins. We test this hypothesis using a fossil-calibrated phylogeny of mitochondrial genomes (mitogenomes) from all extant and recently extinct penguin taxa. Our temporal analysis demonstrates that numerous recent island-endemic penguin taxa diverged following the formation of their islands during the Plio-Pleistocene, including the Galápagos (Galápagos Islands), northern rockhopper (Gough Island), erect-crested (Antipodes Islands), Snares crested (Snares) and royal (Macquarie Island) penguins. Our analysis also reveals two new recently extinct island-endemic penguin taxa from New Zealand's Chatham Islands: Eudyptes warhami sp. nov. and a dwarf subspecies of the yellow-eyed penguin, Megadyptes antipodes richdalei ssp. nov. Eudyptes warhami diverged from the Antipodes Islands erect-crested penguin between 1.1 and 2.5 Ma, shortly after the emergence of the Chatham Islands (∼3 Ma). This new finding of recently evolved taxa on this young archipelago provides further evidence that the radiation of penguins over the last 5 Ma has been linked to island emergence. Mitogenomic analyses of all penguin species, and the discovery of two new extinct penguin taxa, highlight the importance of island formation in the diversification of penguins, as well as the extent to which anthropogenic extinctions have affected island-endemic taxa across the Southern Hemisphere's isolated archipelagos.

Multiscale climate change impacts on plant diversity in the Atacama Desert.

Comprehending ecological dynamics requires not only knowledge of modern communities but also detailed reconstructions of ecosystem history. Ancient DNA (aDNA) metabarcoding allows biodiversity responses to major climatic change to be explored at different spatial and temporal scales. We extracted aDNA preserved in fossil rodent middens to reconstruct late Quaternary vegetation dynamics in the hyperarid Atacama Desert. By comparing our paleo-informed millennial record with contemporary observations of interannual variations in diversity, we show local plant communities behave differentially at different timescales. In the interannual (years to decades) time frame, only annual herbaceous expand and contract their distributional ranges (emerging from persistent seed banks) in response to precipitation, whereas perennials distribution appears to be extraordinarily resilient. In contrast, at longer timescales (thousands of years) many perennial species were displaced up to 1,000 m downslope during pluvial events. Given ongoing and future natural and anthropogenically induced climate change, our results not only provide baselines for vegetation in the Atacama Desert, but also help to inform how these and other high mountain plant communities may respond to fluctuations of climate in the future.

Ancient DNA of crested penguins: Testing for temporal genetic shifts in the world's most diverse penguin clade.

Human impacts have substantially reduced avian biodiversity in many parts of the world, particularly on isolated islands of the Pacific Ocean. The New Zealand archipelago, including its five subantarctic island groups, holds breeding grounds for a third of the world's penguin species, including several representatives of the diverse crested penguin genus Eudyptes. While this species-rich genus has been little studied genetically, recent population estimates indicate that several Eudyptes taxa are experiencing demographic declines. Although crested penguins are currently limited to southern regions of the New Zealand archipelago, prehistoric fossil and archaeological deposits suggest a wider distribution during prehistoric times, with breeding ranges perhaps extending to the North Island. Here, we analyse ancient, historic and modern DNA sequences to explore two hypotheses regarding the recent history of Eudyptes in New Zealand, testing for (1) human-driven extinction of Eudyptes lineages; and (2) reduced genetic diversity in surviving lineages. From 83 prehistoric bone samples, each tentatively identified as 'Eudyptes spp.', we genetically identified six prehistoric penguin taxa from mainland New Zealand, including one previously undescribed genetic lineage. Moreover, our Bayesian coalescent analyses indicated that, while the range of Fiordland crested penguin (E. pachyrhynchus) may have contracted markedly over the last millennium, genetic DNA diversity within this lineage has remained relatively constant. This result contrasts with human-driven biodiversity reductions previously detected in several New Zealand coastal vertebrate taxa.

Subsistence practices, past biodiversity, and anthropogenic impacts revealed by New Zealand-wide ancient DNA survey.

New Zealand's geographic isolation, lack of native terrestrial mammals, and Gondwanan origins make it an ideal location to study evolutionary processes. However, since the archipelago was first settled by humans 750 y ago, its unique biodiversity has been under pressure, and today an estimated 49% of the terrestrial avifauna is extinct. Current efforts to conserve the remaining fauna rely on a better understanding of the composition of past ecosystems, as well as the causes and timing of past extinctions. The exact temporal and spatial dynamics of New Zealand's extinct fauna, however, can be difficult to interpret, as only a small proportion of animals are preserved as morphologically identifiable fossils. Here, we conduct a large-scale genetic survey of subfossil bone assemblages to elucidate the impact of humans on the environment in New Zealand. By genetically identifying more than 5,000 nondiagnostic bone fragments from archaeological and paleontological sites, we reconstruct a rich faunal record of 110 species of birds, fish, reptiles, amphibians, and marine mammals. We report evidence of five whale species rarely reported from New Zealand archaeological middens and characterize extinct lineages of leiopelmatid frog (Leiopelma sp.) and kakapo (Strigops habroptilus) haplotypes lost from the gene pool. Taken together, this molecular audit of New Zealand's subfossil record not only contributes to our understanding of past biodiversity and precontact Maori subsistence practices but also provides a more nuanced snapshot of anthropogenic impacts on native fauna after first human arrival.

Evolutionary factors affecting the cross-species utility of newly developed microsatellite markers in seabirds.

Microsatellite loci are ideal for testing hypotheses relating to genetic segregation at fine spatio-temporal scales. They are also conserved among closely related species, making them potentially useful for clarifying interspecific relationships between recently diverged taxa. However, mutations at primer binding sites may lead to increased nonamplification, or disruptions that may result in decreased polymorphism in nontarget species. Furthermore, high mutation rates and constraints on allele size may also with evolutionary time, promote an increase in convergently evolved allele size classes, biasing measures of interspecific genetic differentiation. Here, we used next-generation sequencing to develop microsatellite markers from a shotgun genome sequence of the sub-Antarctic seabird, the thin-billed prion (Pachyptila belcheri), that we tested for cross-species amplification in other Pachyptila and related sub-Antarctic species. We found that heterozygosity decreased and the proportion of nonamplifying loci increased with phylogenetic distance from the target species. Surprisingly, we found that species trees estimated from interspecific FST provided better approximations of mtDNA relationships among the studied species than those estimated using DC , even though FST was more affected by null alleles. We observed a significantly nonlinear second order polynomial relationship between microsatellite and mtDNA distances. We propose that the loss of linearity with increasing mtDNA distance stems from an increasing proportion of homoplastic allele size classes that are identical in state, but not identical by descent. Therefore, despite high cross-species amplification success and high polymorphism among the closely related Pachyptila species, we caution against the use of microsatellites in phylogenetic inference among distantly related taxa.